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Abstract. There exist many situations where an ordinary differential equation admits a 
movable critical singularity which the test of Kowalevski and Gambier fails to detect. Some 
possible reasons are : existence of negative Fuchs indices, insufficient number of Fuchs indices, 
multiple family, absence of an algebraic leading order. Mainly giving examples, we present the 
methods which answer all these questions. They are all based on the theorem of perturbations 
of Poincare and computerizable. 
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The "usual Painleve test" of Kowalevski and Gambier [H, also called method of pole- 
like expansions, establishes necessary conditions for the Pavrueve property (PP) (defined as 
the absence of movable critical singularities in the general solution of an ordinary differential 
equation). It does so by checking the existence of all converging Laurent series with a finite 
principal part susceptible to represent either the general solution or a particular solution. 

Whenever there exists a movable multivaluedness in the general solution and the test of 
Kowalevski and Gambier fails to detect it, there exists a plethora of other methods able to 
perform this detection. The purpose of these notes is to explain these algorithmic methods not 
in full theoretical detail but on examples selected for their simplicity. 

More details can be found in the lecture notes of a Cargese school ||. 



1 Examples to be settled 

The following three ODEs possess a general solution with movable logarithms undetected 
by the "usual Painleve test" of Kowalevski and Gambier 

u" + Auu' + 2m 3 = 0, (1) 
u"" + 3im" - Au' 2 = 0, (2) 
v!" + uu" - 2u' 2 = 0. (3) 

The first equation (P has a single family of movable singularities u ~ uqx p , X = x ~ x o (one 
should avoid the term branch which induces a confusion with branching i.e. multivaluedness) 

p = -1, E = u (u - l) 2 = 0, indices (-1,0), u = x~\ (4) 

with the puzzling fact that uq should be at the same time equal to 1 according to the equation 
Eq = 0, and arbitrary according to the Fuchs index 0. The Laurent series is here reduced to 
its first term and only depends on one arbitrary constant. The movable logarithms, initially 
found by the a— method (fl2|l §13, p. 221), are exhibited in section 



The second equation (|2J), isolated by Bureau (0 p. 79), possesses two families 

p = -2, u = -60, indices (-3, -2, -1, 20), u = -60x~ 2 + u 20 x 18 + ■■■ (5) 
p = —3, u arbitrary, indices (—1, 0), u — u x~ 3 — 60x~ 2 - (6) 

The first family has enough indices but not enough positive ones, while the second one has 
not enough indices, therefore none of the two families can represent the general solution. The 
movable logarithms are found by two methods, in sections [| and |5|. 

The third equation has no power-law leading behaviour. Chazy ||, |J had to establish a 
special theorem, using divergent series, to exhibit the movable logarithms. The failure appears 
in section |. 

A feature common of the method of Kowalevski when applied to these three ODEs is the 
impossibility to represent the general solution by some Laurent series with a finite principal 
part. Accordingly, the possible presence of multivaluedness in the missing part of the general 
solution cannot be tested by that method. 

The common principle to all the methods performing such a detection is to perturb a 
particular solution into the general solution. Let us first recall the relevant theorems. 



2 Basic theorems 



Boldface letters denote multicomponent quantities. 
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Theorem of perturbations (Poincare, Mecanique celeste [[HI). Consider an ODE of order N, 
of degree one in the highest derivative, depending on a smaTTcomplex parameter e, defined in 
the canonical form 

^ = Kb, u, e], x E C, u e C N , e eC. (7) 

ax 

Let (xo,Uq,0) be a point in C x C N x C and D be a domain containing (xq,uq,0). If K is 
holomorphic in D, 

• there exists a solution u(x,e) satisfying the initial condition u(xq, 0) = Uo, 

• it is unique, 

• it is holomorphic in a domain containing (x , u , 0). 

Proof. See any textbook. Note that K may be independent of e, in which case this is just 
the existence theorem of Cauchy. 

Remark. More practically, the canonical form can also be defined as 

AN 

±£ = K[x,u,u',...,v< lf -% (8) 

Definition. Given a differential equation (DE) 

E(x,u) = (9) 

and a point Uo, the linear DE 

ln] , , E(x,u(°) + Av)-E(x,u(°)) , , 

E' x, u (0) v = lim '- K — '- = 10 

v ' A^O A 

in the unknown v is called the linearized equation in the neighborhood of u^ ) associated to 
the equation E(x, u) = 0. This was introduced by Darboux pj under the name "equation 
auxiliaire". The derivative E' is known under various names : Gateaux derivative, linearized 
map, tangent map, Jacobian matrix, and sometimes Frechet derivative. 
Let us define the formal Taylor expansions 

+oo +oo 

U = en ^ (n \ K = E ^K W . (11) 

n=0 n=0 

The single equation (|7|) is equivalent to the infinite sequence 

rhi(°) 

n = 0: ^— = K(°) = K[x,u(°),0] (12) 
ax 

rlii( n ) 

7i > 1 : ^— = K^ n) = K'\x, u^ 0) , 0]u^ + R^(x, u^, . . . , u (n_1) ). (13) 
ax 

At order zero, the equation is nonlinear. 

At order one, the equation, in the particular important case when K is independent of s, 
is the linearized equation (without rhs, since R/ 1 -* = 0) canonically associated to the nonlinear 
equation. 

At higher orders, this is the same linearized equation with different rhs arising from 
the previously computed terms, and only a particular solution is needed to integrate. 

Theorem //(Poincare 1890, Painleve 1900, Bureau 1939). Take the assumptions of previous 
theorem. If the general solution of (0) is single valued in D except maybe at e — 0, then 



e = is no exception, i.e. the general solution is also single valued there, 
every u( n ) is single valued. 



Proof. See JL2] p. 208. The main difficulty is to prove the convergence of the series. This 
theorem remains valid if one replaces "single valued" (Painleve version) by "periodic" (Poincare 
version) or "free from movable critical points" (Bureau version). 

This feature (one nonlinear equation (]12f ), one linear equation (O) with different rhs) is 
a direct consequence of perturbation theory, it is common to all methods aimed at building 
necessary stability conditions (following Bureau, we call stable an ODE with the PP). The 
equations may be differential like (|i~2j) -( [I3|) , or simply algebraic. Moreover, all the methods 
which we are about to describe (except the one of Painleve) will reduce the differential problems 
to algebraic problems keeping the same feature, and the overall difficulty will be to solve 
one nonlinear algebraic equation, then one linear algebraic equation with a countable number 
(practically, a finite number) of rhs. 

All the methods of the Painleve test are applications of the last theorem : 

1. the method of pole-like expansions of Kowalevski and Gambier ||, 

2. the a— method of Painleve [ 12 ], 



3. the method of Bureau (T|, 

4. the Fuchsian perturbative method 



5. the nonFuchsian perturbative method [11|. 



These methods establish necessary conditions for the Painleve property by building one or 
two perturbed equations from the original unperturbed equation, then by applying the theorem 
II at a point Xq which is movable. This movable point can be either regular (method of Painleve) 
or singular noncritical (all the others), which will require its previous transformation to a regular 
point (by a transformation close to u — > u^ 1 ) for theorem II to apply. One is thus led to the 
equations (|13"1), i.e. to one linear DE with a sequence of rhs. In order to avoid movable critical 
points in the original equation, one requires single valuedness in a neighborhood of xq for : 
the general solution of the linear homogeneous equation, a particular solution of each of the 
successive linear inhomogeneous equations. 



3 Meaning of the negative Fuchs indices 

Basically, they are just the consequence of a resummation of a series. In particular, they 
do not imply the existence of an essential singularity. This is easier to understand if one starts 
from a given general solution rather than from a given ODE whose general solution may not 
be known in closed form. 

The ODE with a meromorphic general solution M 

E = u" + 3uu' + u 3 = 0, u — 1 -, a and b arbitrary, (14) 

x — a x — b 

has two families, 
(Fl) p = —1, u = 1, indices (—1, 1), u = x" 1 + • • •, 
(F2) p = -1, M = 2, indices (-2, -1), u = 2 X ~\ 
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and the negative index —2 must coexist with the meromorphy. The representation of the 
general solution flT4| ) as a Laurent series of x — xq is the sum of two copies of an expansion of 
1/ (x — c) , and there exist two expansions of 1 / (x — c) 

-l 

(x — c)~ l = ^2 (c — Xq)^ 1 ^^ (x — XqY , |c — Xq| < \x — Xq\ (15) 
j=-oo 
+00 

= ~( c ~ ^o) 1- ^^ - x y, \x - x \ <\c- x Q \. (16) 

j=0 

The family (Fl) corresponds to the sum (first expansion with c = a = Xq) plus (second 
expansion with c = b), while the family (F2) corresponds to the sum (first expansion with 
c = a) plus (second expansion with c = b). 

The second family series terminates : u = 2/{x — Xq), and its algorithmic perturbation into 
the doubly infinite Laurent series is handled in Ref. by the method of section (|. 



4 The Fuchsian perturbative method 

It allows to extract the information contained in the negative indices [§], thus building 
infinitely many necessary conditions for the absence of movable critical singularities of the 
logarithmic type ||. 

The perturbation which describes it is close to the identity 

+oo +oo 

x unchanged, u = £ eV n) : E = £ e n B (n) = 0, (17) 

n=0 n=0 

where, like for the a— method, the small parameter e is not in the original equation. 
Then, the single equation (^) is equivalent to the infinite sequence 

n = E (0) =E(x,u (0) ) = (18) 
Vn > 1 E (n) =E'(x, u (0) )u (rt) + R (n) (x, u (0) , . . . , u^) = 0, (19) 

with identically zero. From Theorem II, necessary stability conditions are 

- the general solution u^ ^ of (^) has no movable critical points, 

- the general solution ir 1 ) of (^) has no movable critical points, 

- for every n > 2 there exists a particular solution of (O) without movable critical points. 

Order zero is just the complete equation for the unknown u^, so, in order to get some 
additional information, one must apply Theorem II for a perturbation different from ( U7|) . One 
then uses the method of pole-like expansions at this order zero, only to obtain the leading term 

u(°) ~ u^x p °f a U the acceptable families of movable singularities. 

The precise steps of the algorithm are detailed in Refs. || ||. The algorithm is purely 
algebraic, i.e. one does not perform any integration. The expression of each Taylor coefficient 
u («) i s 

found to be a Laurent series of x — x ~ x o, 

+oo 

u (n) = ^ufV+P, (20) 

j=np 
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in which p is the smallest Fuchs index (therefore lower than or equal to —1), so that the full 
expression u is a "full" Laurent series, i.e. one whose powers extend to both infinities. 

The Fuchsian perturbative method (as well as the nonFuchsian one which will be seen 
section is useful if and only if the zeroth order n = fails to describe the general solution. 
This may happen for two reasons. The most common one is a negative Fuchs index in addition 
to —1 counted once, the second, less common one is a multiplicity higher than one for some 
family, as in the single family of equation ([[]). 



4.1 The simplest constructive example 

The equation ([[]) is the simplest example to understand this method, because 

1. there exists a movable logarithm, as shown by the a— method ([[12] §13, p. 221), 

2. the method of pole-like expansions fails to find it, 

3. the Fuchsian perturbative method finds it after a very short computation, easy to do by 
hand since all series terminate. 

The single family, in the notation of this section, is 

p = -1, = 4 0) («o 0) - I) 2 = 0, indices (-1,0), (21) 
and the order n = yields the one-parameter series 

= x^ 1 (t ne series terminates). (22) 

At order n = 1, the derivative of E at point u = is 

E'(x,uW) = d 2 x + 4 X - 1 d x + 2 X - 2 , (23) 

so that the linearized equation for is of Fuchsian type (by the way, this is why we call its 
indices "Fuchs indices" and not "resonances" , a word which refers to no resonance phenomenon) . 
The computation of the Laurent series ( |2"0| ) for vS 1 ' (with n — l,p — — l,p = —1) is made in 
one computer loop by increasing values oTj and exhibits no logarithms. The result is 

= u^\x~ 2 + u^x 1 ! u -l an d Uq arbitrary. (24) 

The sum u,(°) +eu^' is the beginning of a series which now depends on two arbitrary parameters, 
namely xq and euq (one can indeed set ul\ to zero without loss of generality since it represents 
a perturbation of x ). With the retained value 

= u^x^ 1 , Uq arbitrary, (25) 

the order n = 2 is 

E^ = E'(x,uW)u® +6«( ( V 1 ) 2 +4 M (V 1 )' 

= X- 2 (xV 2 ))" + 2^ 1)2 X - 3 = 0, (26) 

and the Laurent series (|20|) for the particular solution is found not to exist because of a 
movable logarithm 

= -24 1)2 X - 1 (Log X -l). (27) 

The movable logarithmic branch point is therefore detected in a systematic way at order n = 2 
and index % = 0. 

The necessity to perform a perturbation arises from the multiple root of the equation for 
«o°\ responsible for the insufficient number of arbitrary parameters in the zeroth order series 
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4.2 An example needing order seven to conclude 

In the equation @, the first family provides, at zeroth order, only a two-parameter expan- 
sion and, when one checks the existence of the perturbed solution 



u 



n=0 



EM j-2-Zn 
u) 'x J 



3=0 



(28) 



one finds that coefficients i4o , u_\, u_ 2 i u -l can be chosen arbitrarily, and, at order n = 7, one 
finds two violations HI 



Q-l — U 20 U -l 



0,Q 



(7) 
20 



(0) 2 (l) 6 (1) 
U 20 U -3 U -2 



0. 



(29) 



implying the existence of a movable logarithmic branch point. 

Remark [O. The value n = 7 is the root of the linear equation n(z min — p) + (i max — p) = — 1, 
with p = — 2~?min = — 5,i max = 18, linking the pole order p in the Fuchsian case c = 0, the 
smallest and the greatest Fuchs indices. It expresses the condition for the first occurrence of a 
power x 

equation ([19|), r.h.s. created by the nonlinear terms 3mm" — 4m 

As to the second family, it is useless for the Fuchsian perturbative method, because the 
two arbitrary coefficients corresponding to the two Fuchs indices (—1, 0) are already present at 
zeroth order. 



1 1 leading by integration to a logarithm, in the r.h.s. RS n ^ of the linear inhomogeneous 
19]), r.h.s. created by the nonlinear terms 3mm" - 1 



5 The nonFuchsian perturbative method 

Every time the family under study has less Fuchs indices than the differential order N, 
the Fuchsian perturbation method fails to build a representation of the general solution, thus 
possibly missing some stability conditions. Such an example is the second family of the equation 
(fj). The missing solutions of the auxiliary equation (pi]) are then nonFuchsian solutions. 

Although there is no difficulty to algorithmically compute expansions for the nonFuchsian 
solutions, these are of no immediate help, due to their generic divergence. 

There is one situation where some stability conditions can be generated algorithmically. It 



occurs when the two following conditions are met (11 



1. There exists a particular solution u = u^ -* which is known globally, is meromorphic and 
has at least one movable pole at a finite distance denoted x$. 

2. The only singular points of the linearized equation = are x = xo, nonFuchsian, and 
x = 00, Fuchsian. 

Then, the property that a fundamental set of solutions u* 1 ) be locally single valued near 
X = x — x = is equivalent to the same property near x — 00. This is the global nature of u^°) 
which allows the study of the point x = 00, easy to perform with the Fuchsian perturbative 
method. 

An important technical bonus is the lowering of the differential order N of equation E*- 1 ) = 
by the number M of arbitrary parameters c which appear in u^ ^. Indeed, again since u^°) is 
closed form, its partial derivatives d c u^ are closed form and are particular solutions of E^ = 0, 
which allows this lowering of the order. 

At each higher perturbation order n > 2, one similarly builds particular solutions u^ ra ^ as 
expansions near x = 00 and one requires the same properties. 
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Let us just take one example. In section ]4.2| , the fourth order equation (fj) has been proven 
to be unstable after a computation practically intractable without a computer. Let us now 
prove this result without computation at all [II]. For the global two-parameter solution 

= cx~ 3 — 60x~ 2 , (c,Xq) arbitrary, (30) 

the linearized equation 

= E '(x, u®)uW = [&t + 3u^d 2 x - Su^'d* + 3u<°>> (1) = (31) 

has only two singular points x = (nonFuchsian) and x = 00 (Fuchsian), it admits the two 
global single valued solutions d XQ u^ and d c u^°\ i.e. u^ 1 ' = x _4 )X~ 3 - The lowering by M = 2 
units of the order of the linearized equation ([31]) is obtained with 



u 



« = X ~ 4 v : [d 2 x - lQx^d x + Sex' 3 - 60 X ->" = 0, (32) 



and the local study of x — 00 is unnecessary, since one recognizes the confluent hypergeometric 
equation. The two other solutions in global form are 

c^O-.v'l = x^i(24;-3c/x) = X 17/ V 23 ( v /l2^), (33) 
X 17/2 N 23 (Jl2^), (34) 



where the hypergeometric function Fi(24; — 3c/ x) is single valued and possesses an isolated 
essential singularity at x = 0, while the fonction N 2 3 of Neumann is multivalued because of a 
Logx term. 

Remark. The local study of (|3l]) near x — provides the formal expansions for the two 
nonFuchsian solutions 



X 



0, c + : = e+V^^x^il + O(Vx)), (35) 



detecting the presence in ( ^11) of an essential singularity at x = 0, but the generically null 
radius of convergence of the formal series forbids to conclude to the multivaluedness of 
A nonobvious result is the existence, as seen above, of a linear combination of the two formal 
expansions (|35| ) which is single valued. 

An application to the Bianchi IX cosmological model can be found in Ref. [ID|. 



6 Miscellaneous perturbations 

The differential complexity of the a— method explains why it usually succeeds in case of 
failure of all the other methods, which only have an algebraic complexity. 

For equation @, there exists no perturbation satisfying the assumptions of Theorem II 
page f|, there only exist singular perturbations, i. e. which discard the highest derivative. Since 
they however give the correct information, it would be desirable to extend Theorem II in that 
direction. Meanwhile, the reasoning given below cannot be considered as a proof, and one 
should refer to the proof of Chazy. 

Equation @ is handled by the singular perturbation 

+oo +oo 

u = e- 1 J2z n u {n) ,E = e- 2 Y,£ n E {n \ (36) 

n=0 n=0 
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which excludes u'" from the simplified equation and for which all the series happen to terminate, 
which makes the computation quite easy 

E (o) = u (o )u (oy> _ 2u (or = o (37) 

= cx^ 1 , X — x ~ x 07 ( x 0i c ) arbitrary, (38) 
= c( X - 3 (x 2 u (1) r - 6 X - 4 ) = 0, (39) 
Gx-^Logx-l). (40) 



u 



(i) 



7 Conclusion 

One should not be afraid of negative Fuchs indices. Indeed, at the level of the linearized 
equation i.e. of the term n — 1 of a perturbation process, the sign of these Fuchs indices does not 
matter at all. Their presence does not imply the existence of essential singularities, as shown 
on an elementary example, and there do exist algorithmic methods to build no-log conditions 
from them, in the same way as the widely known method of pole-like expansions builds no-log 
conditions from the positive indices. 
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